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Abstract. Using the isospin-dependent quantum molecular dynamical (IQMD) 
model, we systematically study the role of momentum dependent interactions in global 
stopping and analyze the effect of symmetry energy in the presence of momentum 
O ■ dependent interactions. For this, we simulate the reactions by varying the total mass 

2 ■ of the system from 80 to 394 at different beam energies from 30 to 1000 MeV/nucleon 

over central and semi-central geometries. The study is carried in the presence of 
momentum dependent interactions and symmetry energy by taking into account hard 
equation of state. The nuclear stopping is found to be sensitive towards the momentum 
l/^ ' dependent interactions and symmetry energy at low incident energies. The momentum 

CN . dependent interactions are found to weaken the finite size effects in nuclear stopping. 
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1. Introduction 

The heavy-ion collisions at intermediate energies have witnessed several rare phenom- 
ena such as multifragmentation, disappearance of flow, partial(or complete) stopping as 
well as sub threshold particle production [HE]. The recent advances in the radioactive 
nuclear beam (RNB) physics is providing scientific community a unique opportunity to 
investigate the isospin effects in heavy-ion colhsions (HIC's) [3] with respect to the above 
rare phenomena [H |2] . Beside the many existing radioactive beam facilities, many more 
are being constructed or under planning, including the Cooling Storage Ring (CSR) fa- 
cility at HIRFL in China, the Radioactive Ion Beam (RNB) factory at RIKEN in Japan, 
the FAIR/GSI in Germany, SPIRAL2/GANIL in France, and Facility for Rare Isotope 
Beam (FRIB) in the USA [1]. These facilities offer possibility to study the properties 
of nuclear matter or nuclei under the extreme conditions of large isospin asymmetries. 
Though at low incident energies, where fusion and related phenomena are dominant, 
systematic studies over isospin degree of freedom are available [5], no such studies are 
available at intermediate energies. One of the cause could be the much more complex 
dynamics involved at intermediate incident energies. Among various above mentioned 
phenomena nuclear stopping of the colliding matter has gained a lot of interest since 
it gives us possibility to examine the degree of thermalization or equilibration of the 
matter. 

In a recent communication ^ 2\ , nuclear stopping has been explored with reference to 
the isospin degree of freedom. Nuclear stopping in heavy-ion collisions has been studied 
by the means of rapidity distribution [8] or by the asymmetry of nucleonic momentum 
distribution [9]. Bauer and Bertsch [9] reported that, nuclear stopping is determined by 
both the mean field and in-medium NN cross-sections. Different authors suggested that 
the degree of approaching isospin equilibration provides a mean to prove the power of 
nuclear stopping in HIC's [TO]. 

Interestingly, no systematical study is available in the literature on the effect of isospin 
degree of freedom via symmetry energy on nuclear stopping. This becomes much more 
important if one acknowledges that the effect of symmetry energy could be altered in 
the presence of momentum dependent interactions, which has become essential part of 
the interaction for any reasonable dynamical model [HI ^U\. Thus, our aim is at least 
two folds: 

• We plan to understand the role of symmetry energy in the presence of momentum 
dependent interactions in a systematic way. 

and 

• To further examine how mass dependence alter the above findings. It is worth 
mentioning that the study of the mass dependence is very essential for any 
meaningful conclusion. 

This study is done within the framework of isospin-dependent quantum molecular dy- 
namics (IQMD) model, which is discussed in section-2. Section-3 contains the results 
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and summary is presented in section-4. 



2. ISOSPIN-DEPENDENT QUANTUM MOLECULAR DYNAMICS 
(IQMD) MODEL 

The isospin-dependent quantum molecular dynamics (IQMD) [IT] model treats different 
charge states of nucleons, deltas and pions explicitly, as inherited from the VUU model 
[TT] . The IQMD model has been used successfully for the analysis of large number of 
observables from low to relativistic energies. The isospin degree of freedom enters into 
the calculations via symmetry potential, cross-sections and nucleon-nucleon interactions. 
The details about the elastic and inelastic cross-sections for proton-proton and neutron- 
neutron collisions can be found in ref.|ll]. 
In this model, baryons are represented by Gaussian-shaped density distributions 

/,(f,p,t) = 4^ ■ e-^'-""^^ ■ e-(^-'«)^i^. (1) 

Nucleons are initialized in a sphere with radius R = 1.12A^^^ fm, in accordance 
with the liquid drop model. Each nucleon occupies a volume of h^, so that phase 
space is uniformly filled. The initial momenta are randomly chosen between and 
Fermi momentum (p^;'). The nucleons of target and projectile interact via two and 
three-body Skyrme forces, Yukawa potential and momentum-dependent interactions. 
These interactions are similar as used in the molecular dynamical models like quantum 
molecular dynamics (QMD)|8] and relativistic QMD [12]. The isospin degree of freedom 
is treated explicitly by employing a symmetry potential and explicit Coulomb forces 
between protons of colliding target and projectile. This helps in achieving correct 
distribution of protons and neutrons within nucleus. 
The hadrons propagate using Hamilton equations of motion: 
dvi d{ H ) dpi d{ H ) 



with 



dt dpi ' dt dvi 



{H)={T)+{V) 

^2 
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The baryon-baryon potential V^\ in the above relation, reads as: 

V [r — r) = Vskyrme + ^Yukawa + ^Cc 



Skyrme '^ Yukawa '^ Coul 
' mdi ' sym 
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+ t^\n^[t,{p^-p)^ + lW-r) 

+ t,-T^TiS{r-l'-r-). (4) 

Here Zi and Zj denote the charges of i*^ and j*'^ baryon, and T3, T3 are their respective 
T3 components (i.e. 1/2 for protons and -1/2 for neutrons). Meson potential consists of 

Coulomb interaction only. The parameters fj, and ti, , te are adjusted to the real part 

of the nucleonic optical potential. For the density dependence of the nucleon optical 
potential, standard Skyrme-type parameterization is employed. As is evident, we choose 
symmetry energy that depends linearly on the baryon density. 

The binary nucleon-nucleon collisions are included by employing collision term of well 
known VUU-BUU equation. The binary collisions are done stochastically, in a similar 
way as are done in all transport models. During the propagation, two nucleons are 
supposed to suffer a binary collision if the distance between their centroids 



r 



3\ 



O'tot 



< y ^, (^tot = (^{y/s, type), (5) 

"type" denotes the ingoing collision partners (N-N, N-A, N-tt,..). In addition, Pauli 
blocking (of the final state) of baryons is taken into account by checking the phase 
space densities in the final states. The final phase space fractions Pi and P2 which are 
already occupied by other nucleons are determined for each of the scattering baryons. 
The collision is then blocked with probability 

Puock = 1-(1-Pl)(l-P2). (6) 

Delta decays are checked in an analogous fashion with respect to the phase space of the 
resulting nucleons. Recently, several studies have been devoted to pin down the strength 
of the NN cross-section El. 



3. Results and Discussion 

The global stopping in heavy-ion collisions has been studied with the help of many 
different variables. In earher studies, one used to relate the rapidity distribution with 
global stopping. The rapidity distribution can be defined as [131 E 



Y{^) = -/n#l±M) , (7) 

^' 2 E(^)-p,(^)' ^^ 

where E{i) and Pz{i) are, respectively, the total energy and longitudinal momentum of 
i^'^ particle. For a complete stopping, one expects a single Gaussian shape. Obviously, 
narrow Gaussian indicate better thermalization compared to broader Gaussian. 
The second possibility to probe the degree of stopping is the anisotropy ratio (R) |6] : 

''(e.IpiiWI) 

where, summation runs over all nucleons. The transverse and longitudinal momenta 
are p±{i) = JpHi) +P^(0 and p\\{i) = Pz{i)i respectively. Naturally, for a complete 
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stopping, R should be close to unity. 

For the present study, simulations were carried out for the reactions 2f)Ca^^+2oCa^^ , 
28iVi58 +28 Ni^^, iiNb^^ +41 Nb^^, 54^ei3i +54 Xc^^i andjg Au^^'^ +79 Au^^'^ at different 
beam energies ranging between 30 and 1000 MeV/nucleon at central and semi-peripheral 
geometries. The incident energy of 30 MeV/nucleon is the lowest limit for any semi- 
classical model. Below this incident energy, quantum effects as well as Pauli blocking 
need to be redefined. A hard (H) and hard momentum dependent (HM) equation of 
state (EOS) has been employed with symmetry energy Esym = and 32 MeV. The 
corresponding values of the symmetry energy are indicated as subscript {Hq, H32, HMq 
and if M32.). It is worth mentioning that global stopping is insensitive toward the nature 
of static equation of state. Our present attempt is to perform a theoretical study with 
wider variation in the value of symmetry energy and to see how much it can affect the 
heavy-ion dynamics. Though MDI destabilizes the nuclei, a careful analysis is made by 
Puri et al. [15] and found that upto 200 fm/c, emission of the nucleons with momentum 
dependent interactions is quite small. 

Lets start with the aspect of rapidity distribution as a indicator for nuclear stopping. As 
discussed earlier, nuclear stopping is a phenomena which originates from the participant 
zone. To study the nuclear stopping in term of rapidity distribution in Fig{T], we display 
the rapidity distribution of free particles (which originates from the participant zone) 
for different forms of the hard equation of state {Hq, H32, HMq and HM^2)- We see 
that free particles emitted in the central collisions form a narrow Gaussian for the 
heavier system as well as at higher incident energy (say E = 400 MeV/nucleon). It 
is evident from here that nuclear stopping is dominating with increase in size of the 
system as well as at higher incident energy (say E = 400 MeV/nucleon). With increase 
in incident energy, this behavior is not expected to be universal. It is clear from the 
Ref.[T6] that maximum stopping is observed around 400 MeV/nucleon. This study is 
done in detail in term of anisotropy ratio R in the Figs. [2] and [3] of letter and found 
to be in supportive nature with Ref. [16]. On the other hand, nuclear stopping in term 
of rapidity distribution of protons is found to be weakly sensitive towards symmetry 
energy and momentum dependent interactions. This is due to the reason that there are 
not only free particles which originates from the participant zone. The other candidate 
which originates from this zone are the light charged particles (LCP's) (2 < A < 4). 
In our recent communication [17], these LCP's are shown to be more sensitive towards 
symmetry energy due to pairing nature and also good indicator for nuclear stopping. 
For more meaningful, we will check the sensitivity of symmetry energy and momentum 
dependent interactions on nuclear stopping in term of anisotropy ratio R, which is 
defined in detail in the above paragraph. In Fig. [21 we display the time evolution 
of the anisotropy ratio R for the central collisions of 2oCa^^ + 2oCa^^ (left panel) 
and ■jgAu^^'' + ■jgAv}'^'^ (right panel). The incident energies of 30, 50, 400 and 1000 
MeV/nucleon are employed using the four different possibilities of hard equation of state 
(ifo) -f^32, HMq and HM^2)- Interestingly, the anisotropy ratio /?, though, is insensitive 
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Figure 1. The rapidity distribution ^ as a function of reduced rapidity for free 
nucleons for the reactions of 2oC'a'"^ +20 Ca^^ , and rgAu^^^ +79 Au^'^'^ for four different 
possibilities of hard equation of state. The top and bottom panels are at E = 50 and 
400 MeV/nucleon. 



towards the symmetry energy, shows appreciable effect for the momentum dependent 
interactions. At high enough incident energy, both effects (of momentum dependent 
interactions as well of symmetry energy) wash away. For further test, we simulated two 
reactions (i) 2oCa*^ +20 Ca^^ and kept the same N/Z ratio by taking 100-^^"^° +100 X^°° 
reaction. In other case (ii) we took neutron rich reaction ygAti^^'^'+ygA-u^^^ and simulated 
the reactions of i%X^^ +16 X^^ by keeping N/Z ratio again same. In both the cases, 
the above trend holds good. Therefore, indicating that the above behavior is universal. 
Even rapidity cuts to mid-rapidity region do not yield different results. Below the beam 
energy of 50 MeV/nucleon, collision dynamics is governed by the mean field. Therefore, 
interactions involving isospin particles like nn, np, pp dominate the outcome and hence 
symmetry effects are visible. 
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Figure 2. The time evolution of anisotropy ratio R for the reactions of 20^0^° +20 
Ca^^ , and rgAu^^^ +79 Au^'^'^ . The panels from top to bottom represent the scenario 
at beam energies 30, 50, 400 and 1000 MeV/nucleon, respectively. 



As discussed earlier, R approaches to 1 at low incident energies, indicating the isotropic 
nucleon momentum distribution of the whole composite system. The behavior at 30 
MeV/nuceon is little different due to the fact that binary collisions do not play any role 
and mean field will take larger time to thermalize the colliding nuclei. As beam energy 
increases above the certain energy, R starts decreasing from 1 towards the lower values, 
indicating partial transparency. This value of the beam energy, above which R starts 
decreasing, depends on the size of the system. In our observations for the reaction of 



79 



Au 



197 



79 



Au \ it is close to 400 MeV/nucleon. This finding is similar to the one 
reported by W. Reisdorf et al. [16]. The value of -R > 1, can be explained by the 
preponderance of momentum perpendicular to beam direction [18]. This is true for all 
equations of state. It is also seen that relaxation time decreases with the increase in the 
beam energy, while, increases with the increase in the mass of the colliding system. It 
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Figure 3. The final state anisotropy ratio i? as a function of composite mass of the 
system Atot for different possibilities of hard equation of state (discussed in text). The 
left and right panels are at central and semi-peripheral geometries. The panels labeled 
with a, b, c and d are at E= 30, 50, 400 and 1000 MeV/nucleon, respectively. All the 
curves are fitted with power law 



shows that higher beam energies and hghter systems lead to more violent NN collisions 
and faster dissipation. This is consistent with the isospin equilibrium process as shown 
by Li et.al. [3]. 

It will be further interesting to see whether the above findings have mass dependence or 
not. This is particular important since the role of the momentum dependent interactions 
and the symmetry energy depends on the size of the system. For this, we display in 
Fig. [31 the anisotropy ratio i? as a function of the composite mass of the system 
{Atot = At + Ap) at different beam energies ranging from 30 to 1000 MeV/nucleon. 
The left panel represents the results at central geometry, while, right panel is at semi- 
peripheral geometry. 
Our findings are: 

• The anisotropy ratio R increases with the composite mass of the system. This is 
true for all incident energies and impact parameters. This dependence becomes 
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weaker as one moves from central to semi-peripheral geometry. It is due to the 
fact that nuclear stopping is governed by the participant zone only. This is further 
supported by the fact that at higher incident energies e.g. E = 1000 MeV/nucleon, 
the stopping is almost independent of the composite mass of the system at semi- 
peripheral geometry and almost 50% decrease is observed in the nuclear stopping 
as compared to central geometry. 

• The effect of the symmetry energy is visible below the Fermi energy. Same 
conclusion was also reported in Fig. O The effect of the symmetry energy 
diminishes in the absence of momentum dependent interactions. Moreover, this 
effect weakens at semi-peripheral geometries. 

and 

• The importance of momentum dependent interactions is also visible in the nuclear 
stopping. This effect decreases as one moves from the low to higher incident energy 
and from central to peripheral geometry. At the higher incident energies e.g. at E 
= 1000 MeV/nucleon, the anisotropy ratio is independent of the equation of state 
and symmetry energy. It is worth mentioning that the inclusion of momentum 
dependent interactions is found to suppress the binary collisions and as a result is 
found to affect the sub threshold particle production as well as disappearance of 
collective flow [H [H] • 



4. Summary 

Using the IQMD model, we have studied the nuclear stopping for isospin effects. For 
this, we simulated the reactions of 2oCa'^° +20 Ca'^^, 2&Ni^^ +28 Ni^^, ^iNlP^ +41 NlP^, 



1 

54X6^^^ +54 Xe^^^ and ygA-u^^^ +79 A-u^^^ in the presence of momentum dependent in- 
teractions and symmetry energy. The role of symmetry energy at low incident energy 
gets enhanced in the presence of momentum dependent interactions. Further, we can 
conclude that maximum stopping is obtained for the heavier systems at low incident en- 
ergies in central collisions in the absence of momentum dependent interactions implying 
that momentum dependent interaction suppresses the nuclear stopping. 
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